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Abstract 

Models of the mai n geomagnetic field are generally represented by a 
scalair potential 7 expanded in a finite number of spherical haunaonics. 
la the last decade, such models have been derived mainly by a recursive 
iteration method from the field, magnitude F observed by satellites in 
low-altitude polar orbits. Very accurate observations of F were used, 
but indications exist that the accuracy ol models derived from them is 
considerably lower. One problem is that F does not always characterize 
7 uniquely: Backus has derived a class of counterexamples in which two 
different choices of 7 correspond to the same F • It is not clear 
whether such ambiguity can be encountered in deriving 7 from F in 
geomagnetic surveys, but there exists a connection, due to the fact that 
the counterexamples of Backus are related to the dipole field, while the 
geomagnetic field is dominated by its dipole component. If the models are 
recovered with a finite error (l.e. they cannot completely fit the data 
and consequently have a small spurious component) this connection allows 
the error in certain sequences of heumbnlc terms in 7 to be enhauiced 
without unduly large effects ^n the fit of F to the model. Computer 
simulations have demonstrated uhis effect, piroducing as a result 
models which fit the data of F quite closely but yield much poorer 
fits to the direction of the magnetic vector. Possible remedies are 
discusc:.d. An appendix also discusses a particulau: claiss of fields - 
related to the counterexamples of Backus - for which it can happen that 
the recursive iteration deriving 7 from F does not converge to the 
correct solution. 
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INTRODUCTION 


In order to derive a general vector field in space, three independent 
scalars must be observed at each poini - e.g. 5 orthogonal components, 
or the field's magnitude and tvo angles defining its direction. For the 
main geomagnetic field B - i.e. the part vhich originates deen inside 
the earth - such an observation is redundant, since that field is comple- 
tely defined by a scalar potential 7 

B = - 77 (1) 

The information needed to derive B "*»y thus be given by a single 
scalar function of position. In fact, there exists the additional con- 
straint that 7 must be harmonic ; this may be viewed as reducing the 
"information content" of 7 to that of a two-dimensional scalar function, 
since only 2 indices appear in the spherical harmonic expansion of 7 , 
traditionally written as 

m=n 

7 = a JF*, (a/r)^*^ ^ { ®n ^2) 

n*l m=0 

( a is the esurth’s radius). In satellite surveys it is relatively easy to 
measure accurately one of the three scalars associated with the field, 
namely the magnitude |B| , commonly denoted in geomagnetic literature 

by the capital letter F . Accurate observation of F requires neither 
information on the attitude of the observing spacecraft, not observation 
of the flexing of the boom carrying the magnetometer ; such information _is 
required for accurate determination of the field’s direction, which con- 
sequently is much harder to perform. 

Furthermore, satellites conducting surveys of the field genereilly do not 
vaiy their altitude by great amounts, so that their observations of F tend 
to bo distributed more or less over a two-dimensional surface - that of 
a sphere concentric with the earth. The preceding qualitative arguments do 
however suggest that the amount of information thus gathered may suffice to 
determine the field. 
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This has led to attempts to trade one scalar for the other, i.e. to 
derive y from observations of F (of course, observed data must first 
be "cleaned" of fields of external origin, but this irapoitant point 
is not considered here). Zmuda ^ ig^S j appears to have been the first 
to propose this approach and its first major use occured in the analyffis 
of Vanguard 3 data by Cain et ^ I962j . The derivation of models of 
the main field from satellite data has generally been based on F 
^ Heppner . 1963 ; Cain et al. , 196? J Cain , 1971 J ^ particular the 
models derived from observeticn.' of OGO 2,4 and 6 ^ Cain and Sweeney . 

1970 j • The method used in such work is described further on . Because 
P ;.s eeisily observed down to intensities of 1 7 ( = 10 gauss ) this 

hat become ths method of choice for mapping the geomagnetic field and models 
which fit observed F within an r.m-s. deviation of 4-77 have been 
derived by it. 

However, two problems arose - one experimental, tte other theoretical. 
Experimentally, the few observations in which the vector direction of the 
field could be compared to the model ixidicated considerable deviations 
r Cohen. mk ; Woodman . isa] . Al so, the results of airborne surveys 
indicate (J. Cain, private communication) that field components often 
deviate by 100-200 7 from the values prescribed by models. 

The theoretical difficulty was pointed out by G. Backus f 1968 , 1970 1 
who investigated the question of whether F observed on a sphere uniquely 
defines the scalar potential 7 from which it can be derived. If the 
expansion (2) contains a finite number of terms it may be shown f Backus . 

that F and 7 are uniquely related. However, if an infinite 
number of terms is allowed, that assertion turns out to be false and 
Backus [M9] furnished a class of counterexamples in which two quite 
distinct infinite expansions, similar tc (2), yield identied forms of 
F on the unit sphere r = 1 . 

This discovery has cast seme doubt on the assertion that the geesnag- 
netic scalar potential 7 is uniquely derivable from its corresponding 
P # So far, this question has not been resolved. 
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What the present vork shows is that the experimental discrrjpancies 
between vector field observations sz'.d the predictions of the aodel aay 
have a mathematical origin, related to the v rk of Backus* Specifically, 
it will be shown that when perfect recovery f 7 is not possible - 
due to finite observaticmal erre s of varic origins - and each of the 
coefficients of (2) contains a certain error, then this error is enhanced 
for certain seg^uences of tenas related to those derived by Backus. The 
basic reeison for the connection is that the counterexas^les constructed 
by Backus are related to the dipole field, while the geomagnetic field 
is known to be dominated by its dipole component. Con^uter simulations of 
the effect and possible remedies will also be described. 


TLE COUNTEREXAMPLE OF BACKUS 

Backus [iSffi] sou^it to determine whether it was possible, in 
three dimensional space, to fizd two harmonic functions . 

vanishing at infinity and differing hy more than their algebraic signs, 
such that on the unit sphere r »= 1 

( 5 ) 


Be found it convenient to reformulate the problem in terns of the sum 
u and the difference v of ttese functions, which on r * 1 have to 
satisfy 


7u • =0 ( 4 ) 

Choosing u to be the dipole potential 

u ■ 7^ “ cos9/r® (5) 

Backus assuned v to be an infinite series of form (?) and he 
then examined the coefficients to see whether they could be nade to 
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satisfy (!{■) . Not only did this turn out to he possible but a whole class 
of such solutions was found to exist. The general form of these solutions 
is best expressed if one interchanges the sunmations of equation (2) 
and writes 



n=ni 


where for sake of definiteness the leading term in each of the summations 
over n is set equeu. to unity 


= = 1 (7) 

m 

Then, provided the series in the inner summations are suitably chosen 

(there exist two such summations for each value of m , and they are 

actually the same - i.e. bT » G^) , equation (4) is satisfiec* for any 

eurbitrary choice of the coefficients a and b . The series over n 

mm 

are derived by recursion from their leading terms, and since these terras 
are fixed by (7), the series are unique. They are described in more 
detail in Appendix I and contain nonzero terms only for values of n 
and m which sum up to an even number. 

Strictly speaking, (6) satisfies (4) only if an infinite number of 
terms is used. We shall denote this limit of v by a superscript 00 
and write 

V7^ * Vv** » 0 (8) 

If the number of terms in v is finite, (8) no longer holds; how- 
ever, the scalar product contributed by such a truncated v (denoted 
by a superscript f for "finite") may be small if v 




is approached 
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sufficiently closely. This may he : j:r' .lljallv vritton as 

Vi'd • = 0(e) (9) 

vhere 6 1 is a small parameter* In vhat follovs ve indicate for 

brevity all small quantities by £, ; at least 3 distinct quantities cf 
this sort appear, but devoting to of them a snarate notation morel;)" 

encumbers the expressions and yields no nev results# The magnitude of 
the right hand side of (9) depends on the exact form of - in 

particular, on the number k of terms in the shortest summations over 
n which appear in it . If k is sufficiently large, however, the 
expression (9) niay be made arbitrarily small. 

MODE! DERIVATIOH FROM SCALAR DATA 


In the derivations of models of the main geomagnetic field from a set of 

observations of F “ |V7| > the value of 7 is determined by iterated 

linearization, starting from an initial approximation 7^ of the scalar 

potential. Let £7 be the correction to be added to 7^ in order 

1 ^ 

to obtain the true potential 7 

7 = 7 «+ £> ( 10 ) 

o 1 

Then 

ZV7^’V7q " 2 ^ ^ ■ 2 ^ 

2 

If the term of order £ is neglected and 7 is expanded in a 
finite series of form (2) with N unknown coefficients, then (ll) 
represents one linear equation with N unknowns for each point at 
which an obser”'’ation of F was made. 

Typically the number M of such observations greatly exceeds N . If, 
for Instance, the series (2) for 7 is truncated past n = , then 

^1 
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N 


n 

3ax 


1 )‘ 


- 1 


(12) 


Typically, ^ > ■while the number of equations (one 

for each observation of F) may be of the order of 10,000 . The best- 
fitting set of coefficients for y is then derived by least squares 

Z. 

fitting, after which one replaces by. ( + € 7 ) and repeats 

the process, seeking a further correction term y' • This procedure is 
continued until the r.m.s. difference between observed values of F 
and those obtained from the model no longer shows any improvement. 
Various refinements may be included, e.g. dividing each of the equations 
(ll) by |7 7Qj order to avoid giving greater weight to observations 
in polar regions, where F is larger. 


Cases are known in which this procedure does not lead to the correct 
solution (Appendix II), but in analysis of simulated data from models 
representing the geomagnetic field this method tends to converge rapidly 
- to the “correct" solution, if the derived series is capable of ex- 
pressing it, or to a small range of final parameters corresponding to 
u '*best-fitting series" if the data contain unresolved higher harmo.iics 
or other sources of small error. 


ERROR ENHANCEMENT 


The basic equations used in the iteration are 

E77o» V7^ = f ^ 

with one such equation being contributed by each observation point. 

These equations will lack a unique solution if the set of equations 

V 7q • V 7^ = 0 (ll*) 

evaluated for the same points, does have a solution. In that case, given 
one solution of (l5) , new solutions may be obtained by adding to it 
arbitr/tiy multiples of solutions of (14) . 
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If all o^)cerv'ations vcre conduct od at the uane radial dict?.n'0, if 

y were the dipole field and if 7 contained infinite n^^ler of 
^ a 

terms ^ then hy (8) there vculd have existed a solution to (l4) and 

(15) vould not have defined 7 unique* /. !Ione cf these conditions 

1 

holds exactly in the analysis of the neoma^ncti? field (in fact. 7^ 

' vj 

changes from one iteration to the next), -but they are all approxlr.ately 
valid and this, in the presence of finite erre may cause problems. 

In this section the finite spread in the radial distance r will 
be neglected (its effects are considered separately later onV Ir..t 7^ 
now be separated into a dipole part and a non-dipole part 


>■0 * >'d + 6 >'nd 


(15) 


where is yet another small quantity# 

If an exact solution for equations (15) does not exist - due to 

truncation, experimental error, unresolved contributions of external 

field sources etc. then the **best fitting” result for 7 will consist 

1 

of the sum of the ”true solution” 7. and an ”error” 7 , each of which 

t e 

will be expressed by a finite series similar to the one in equation (2) • 
One can then write 


vhere - will signify %est fitting "by the least squares method" (for 
our particular set of observations). 

IJow a finite series y of form (2) can in turn be resolved into 

e 

two parts; a sum 7^^ of "Backus terms" similar to those in (6), plus 

a sum 7 . of independent harmonics. To perfon.; such a separation one 
ei 

first isolates all terms with n = m and Includes them In 7^^ , while 
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the remaining terms are assigned to 7 . Each of the terms in 7 . is 

el eh 

now viewed as the leading c.erm of a Backus-type series of the type listed 

in equation (6) an' in Appendix I (including the coefficient a or 

m 

multiplying such a series), and appropriate higher-order terms are 

added to it to form the rest of its series, up to the highest n and 

m allowed by the finite expansion of 7 , ; to maintain the balance, 

et> 

the same terms that are added to 7 , are subtracted from 7 . . 

eb ^ei 

Finally, those expansions in 7^^ deemed to have (after truncation) 
too few terms fcr (9) to be valid (with the right hand side small enough) 
are transferred from 7^^ to 7^^ . This will certainly include all 
t-hose Backus series which contain only one term - and perhaps also those 
containing two . Equation (16) now may be written 

’’’el' " 2 


The cor-^ ributions of the error terms to the left hand side are mainly 

duo to scalar products with the dipole term, which overshadows the 

non-dipole term. If alT coefficients involved in 7^^ and 7^^ were 

equal, then the contribution to the scalar product associated with a 

coefficient of 7^^ will on the average be one order in £ smaller 

than the contribution associated with one of the c 'efficients of 7 . , 

el 

since oy equat?Lon (9) tne former conti but Ion gains ein extra order in £ . 

Alternatively, if the least squares optimization adjusts all L^ources 
of error so that t-hey contribute approximatel y equally , then the coeffici- 
ents of 7^^ will be larger by a factor of the order of £ than those 

of 7 , leading to an enhanced error in these terms. 

e3 

A different view of the same effect is gained by writing the basic 
relation prior bo linearization 








( 18 ) 
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When the left hand side is expanded in orders of t , the tern 

2 

is simply (^ 7 ^) and all 0 (€ ) terms of the scalar potential contri- 
bute to the term of order t , except for 7 ^ , r’hich due to ( 9 ) only 
contributes at the £ level . Thus j, ~ is inher“Mt.3v less sensitive 
to fluctuations of 7 ^^ , allowing il s errors to be relatively large 

without excessive effects on the fit of > Intuitively, the Backu. 
series may be viewed as modes at which a near-dipole field, derived 
the manner described here, prefers to fluctuate - Just as a bridge or a 
beam tend to oscillate in certain modes related to their geometrical 
properties. 

SIMULATION 


A test of the preceding argument was performed as follows, A 120- term 
expansion of 7 (n = 1C was used to derive the corresponding values 

1H3X 

of F c,t iMi-C points on the unit sphere. A program analyzing these 
simulated data by the previously described method was then applied, but 
the derived model contained onl^ 6^, 80 or 99 terms = 7 , 8 or 9 ). 

UIcUC 

Up to 10 iterations were performed, a3.though the model generally showed 
little variation after the second iteration; the results cited here 
therefore refer to the fit after 2 iterations. As a check a differert 
Input for 7 was used and similar effects were observed, although they 
acre not Included in the results listed here. 

Fi’om the models, tlie predicted intensity was derived for the 

14A0 input point and at each point the deviaiion 

5 F = I F - F^ ( 

was computed. The quantity Sf/F was then derived, its average was 
also found and a i-epresentat ive table of its values at every 5 rd point 
was printed out, from which the worst (i.e. largest) value was picked 
visually. In addition, the mean angle (in radians) between the model 
field and the true field was lerivedt it is of the same order as the 
relative discrepancy ^3^/ in one of the components of B . As before, 
the worst case was also picked out visually from a representative tabic. 
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The results are given in Table 1 • The relative accuracy of the direction 
is more than 10 tim^s vorse than that of the magnitude, and the vorst. fit 
in direction has an error about 10 tines larger still* This is an effect 
slightly* larger than vhat is founc. vith models of the geomagnetj.c 
field. 

Inspection of the coefficients reveaJLed that the error in the derived 
models indeed contained relatively large contributions from the Backus 
sequences* For this purpose the differences ^ ^ ^^tveen 

the harmonic coefficients of the model potential and those of the poten- 
tial used as input were tabulated* The coefficients sequences 


g 


m 

‘m 


m 

2 ^ 


m 


and 


m 


m 


iti 


2 ' eH-4 


corresponding to terms which par^ icipate in a single Backus series aliri 
invariably showed the characteristic eOLtemation of signs between conseca- 
tive ^‘^rmr which is typical of any such series* Furthermore, when the 
ratio R of any such term to the leading term of its series is compared 
to the corresponding ration for the Backus series (Table 4) , the 

two ratios are roughly similar* 

Re'^ults from the test for some of these sequences, as veil as for 

some terns not belonging to any sequence, are shown in Table 2 . As can 

be seen, some of the sequences grow to relatively large amplitudes 

( e*g* the series for 99 term recovery, the g^ series for 80 term 
2 3 

recovery) , others ^e.g* the one headed by in the 99 case) 

2 

fail to develop. 
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POSSIBLE REMEDIES 


The coanterexaraplet: discovered hy Backus satisfy equation (4) on.ly on 
a single spherical surface. If observations are distributed over a finite 
region in three dimensions, the series (6) no longer satisfy (li) over the 
entire region of observation and it is thus possible, in principle, to 
overcome the problem described earlier. In fac^ Backus proved [ Backus . 
Mi] that equation ,(l4) cannot be satisfied over a spherical shell of 
any finite xhlckness. 


In practical observations such a three dimensionality alvays exists 
due to the finite eccentricity of the orbit of the observing sateOJ-ite. 
For those satellites of the OGO series vhich were used for geomagnetic 
surveys this eccentricity was rather small: the altitude range (above sea 
level) was about 400 to I5OO km for OGO 2 , 400 to 9OO km for OGO 4 and 
400 to 1100 km for OGO 6 . Unfor’v unate ly, it appears that in the presence 
of finite error such a variation in altitude is only able to reduce some- 
what the effects described ixi this work and will not eliminate them 
altogether. 

Suppose observations are confined to a spherical shell with 


^r ^rr^^Xr 


If Sr/r = 0(6 ), then equation (9) - which is the one relevant 

here - may still bo used, for by Taylor expansion 


u(r, e,'e) = u(r^, e,'») + o(fe) 

v(r, 0 ,'T) = v(r^, e,'f) + 0 ( £ ) 


(19) 


If u and V satisfy (9) on the sphere r =» r^ then "the above 

equations show that they satisfy a r milar equation over the entire 
spherical shell, although the "smallness” of the right Inrid side is 
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somewhat impaired. Th’.:? '-quations (I 7 } and (l^ rtiU. hold and errjr 
enhancement still exists, although net as '■trcxigly as; In the case of 
obser/ations in an infinitely thin spherical shell. 

As a test, the recovery of the 120- term model, used in 

tehles 1 and 2 , was repeated with data points oscillating over a finite 
range in radial distance, comi>aLrable to that of the OGO orbits. A distinct 
improvement in the angular resolution was obtained (Table 3 ) but there 
still exists a pronounced discrepancy between Jf/F and X>. ( the 

values of ^f/f deteriorate somewhat with increasing Sr, probably 
because the "error" is an unresolved set of n = 10 terms, the contribution 
of which increases rapidly for those points at which r is diminished). 

If an analysis similar to that of table 2 is performed, it is found 
that the Backus series are still enhanced, although the ratios R depeurt 
further from those in table it- . 


More generally, the addition of some vector data to the data set cein 

in general help resolve the problem. In such cases it might be useful to 

separate the model potential y derived from P into two parts 7 , 

m jno 

and 7 . , in a similar way to the resolution of 7 in ( 16 ) into its 
HU. e 

two parts r , and 7 . in (l?)* One may then etssune 7 , to be accu- 
eo el ml 

rately derived frern F end use the added data for deriving 7 ^.^ alone. 


Suppose for example that one is given as added vector data an accu- 
rate location of the dip equator on the surface r = a « On the dip 
equator '^y./'^T = 0 and therefore, by ( 6 ) (noting the equality of the two 
summations when ( 7 ) is assumed) 




r=a 


E 

m=l 


(a cos mY + *b sin mY ) / . (n4 l) G 


m 


m 


E 




r 

n 




( 20 ) 


Here the left hand side is known, the coefficients are given in 

table 4 and the unknowns to be derived are the factors a_ and b . If 

zn nz 

the variation of 0 is neglected (assuming it to be neeur 31 / 2 ) the solution 
of (20) reduces to the expansion of a given function of V In a Fourier 
series. In practice, with 0 kept as variable, the sln 5 >lest way of deriving 
a^jj and b^j, would be by least squares fitting of ( 20 ). 



APPENDIX I : 


TEE BACKUS COUNTEREXAMPLE 


The published fora of the counterexajaple of Backus Backus . M2] 

normalizes the Legendre fimctions in the rnanner preferred by 

matheniaticians and expresses the dependence on the longitude angle Y 
in terms of complex exponentials. For computer progi-ams it is much 
more ccnvcrd-cnt to adop* Gaussian normalization, vhich avoids the 
need to calculate square roots when deriving either sidierical harmonics 
or recursion coefficients, Ic is also more ccgiventional in geomagnetic 
usage to express the dependence on Y by means of cos m*? and sin a'f 
(as in equation 2) instead of exp ( * i m ^ ). The purpose of this 
appendix is to list some properties of the Backus counterexamples when 
these are e3q>ressed in Gaussian normalization with trigonometric 
coefficients In H* , 

The general form of the coxmterexample is given in equation (6), 
and if condition ( 7 ) holds, then the two series there are identical - 

i,e. ^ , The basic recursion relation for either of these series 

n n 

is, with Gaussian normalization. 


- (m»2kH2m-» 4k »l)(2m* 4k - l) (A-I-l) 

m-*>2k ” 12k (m»2k + l)(m + k) '*m + 2k ~ 2 


A similsur recursion holds for rf** ; if ( 7 ) '’s assumed, the results 

of this recursion are shown in Tabls 4 , Note that in all terms involved 
here the sum of Indices is even; terms for which this sum is odd vanish 
identically. 

Table 4 shows that for any given ra , the magnitude of the coeffici- 
ents increases with increasing lower index a , Kie resulting series 
nevertheles.. do converge (cs was shown by Backus), because the normali- 
zation factors entering also contribute to the convergence. For 

Gaussian normalization the for any given m are most conveniently 

generated by a recursion relation starting from 
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p®( e ) = sin^e 

in 

.(0) = COS0 sin® 9 

HI T J. 

The relation is 

A, (6) = cos9 I^(e) 

n+ 1 n n n-i 

where 

= (n m)(n - m)/r(2n-f l)(2n - 1)3 

a ► 

The derivatives are obtained by recursions found froa the above 
by differentiation, starting with 

dj^d® = o cos© sin*"^6 


(A- 1-2) 


(A- 1-5) 


(A- 1-4) 


(A- 1-5) 

ones 


(A-I-6) 



APPENDIX II : CASE IN WHICH ITERATION FAILS 


it vas noted that in practice the iteration ( 15 ) tends to converge 
rapidly to an "optimal model" of the field ; the purpose of this 
appendix is to point out a case in which this is not true. The exemple 
given is closely related to the work of Backus ^ 1970 j and is probably 
rather atypical, but it does show that ( 15 ) is not universally useful. 

Let 4"* and 9 ** be two potentials of the Backus type^ 
satisfying on r = 1 the condition 

^ (A-H-1) 

Let and be tznmcated versions of these series, with all 

'1 '2 

terms having upper indices exceeding n omitted, and let 

(A-II-2) 

By the theorem of Backus proving th& uniqueness of the relation 

between 7 and F for potentials with a finite number of harmonic 

terms ^ Backus. 1968 ^ , uniquely characterizes ^ ^ and 

uniquely characterizes • However, as n-^«o , both F*' and F° 

'2 12 

tend to F and therefore, for large values of n , the difference 
between the two becomes rather small. 
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This suggests that If the input to the recursion (ij) consisted of 
or F*^ , vith n large enough, the recursive iteration might 

X 2 

n 

experience difficulty in telling the tvo apart - e,g., given F , it 

1 

might start converging to either or , depending on the 

initial choice of 7 q • 

This is confirmed by computer simulations vith F generated by 

48-term expansions and the least-squares iteration capable of resolving 

an equal number of coefficients. If the program deriving the model of 

7 is presented by such an F^ , one of tvo things is found to occur. 

If the initial approximation resembles more than it resembles 

4 ^ ) the recursive iteration vlU rapidly converge to the correct 
' 2 

coefficients. If, however, the initial approximation is nearer to 4” , 

' 2 

the coefficients vlU tend to evolve tovards those of 4 ^ • The itera- 

^ 2 

tion then never converges - the search being ccxiducted in the vrong 

range of parameters - and Instead the values of the output coefficients 

wander vithin a finite range around those of 4 ^ • 

' 2 

For a recursion deriving 48 term expansions, each resulting 7 
m; •/ be vieved as defining a point in a 48-dlmension6d space in 

which ea^fc coefficient corresponds to one coordinate, and (15) may be 
v:’eweu as prescribing a mapping of each point in another 

point there. If the field magnitude used as input is ( n 4 6 ) » 

the mapping will have a fixed point which maps onto Itself, 

and in some region surrounding this point the mapping vlU converge to 
• However, this convergence evidently does not hold for all of 
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, since for some initial points the process leads to the neighbor- 
hood of the "false solution" A ^ . The region frop vhlch the iteration 

• 2 

converges to thus has a boundary, vhich vill represent a limit 

cycle - i.e. if the initial i>oint is located on that boundary, the 

mapping causes the resulting set of coefficients to vander along the 

boundary without either converging to or approaching . 

T 2 


nonzero 


In one simulation the^^coefficients of were 


0 


10 


g^ 2 


g^ = - 2.8 

4 


^ = 3.375 

6 


(for the corresponding , the signs of the m = 2 terms are 

reversed). The coefficients of the initial field were chosen as 


S-. 


10 


e 


g^ - - 2;» g® « 2.8A 

2 4 

are 

For A = 1 the last two coefficients of the initial field^equaJ. 

coefficients of the "false solution” 4 ^ “id the recursion in that 

' 2 

case heads for that solution (the spurious term g^ is whittled down 
in the process). On the other hand, if ^ = 0 the recursion heeids 
for the "true” solution transition, corresponding to a limit 

cycle, occurs near X - 0.102 and when that value was used, as many 
as 6 iterations were completed before a clear trend became evident. 
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CAETIONS TO TABLES 

Table 1 Deviations of a model geomagnetic field^derived by means of 
P ^frora the "true” field it is supposed to represent, in a com- 
puter simulation. The input field is given by a scalar potential 

vith 120 coefficients (n = lO), the niariber of points is l440 

max 

and 2 iterations are used, except for one case where the 
results of 5 iterations are also shown for conqparison. All 
results are in units of 10 • 

Table 2 Some of the differences Ag “ g (model) - g(true), where g 
stands for one of the harmonic coefficients ^ or h™ in a 
simulated recovery of 7 fr(xn F f?^escribed in Table 1 • 

Here denotes the ratio between a term in a Backus sequence 
having the type and indices indicated on the table and the 
leading term of its sequence (these ratios are lisl id in Table If), 
while R is the ratio between correspcmding values of Ag • 

Table 5 Results similar to those of Table 1 , with a similar input 

model and 99-term recovery in 5 iterations, for cases in which 
data points are spread out over a spherical shell bounded 

Table 4 Ratios derived from the recursion relation A-I-1 • 

xa 2K m . 



No. of terms 
In recovery 

65 

8c 

99 

2 Iter. 

5 iter. 1 

I 

Mean Sf/F 

26 

20 



Worst Sf/F 

2l|-9 1 

207 

35.0 

51.4 

Mean SMradianrs) 

500 ! 

222 

102.8 

106.4 

Worst S/\(radians) 

1 _ . . _ 1 

2497 

i 

1 

2159 

1296 

1 

1507 


TaMe 1 







65 terms 



+> >> A « 


oottOA: tio^A toboA 


Table 





Mean £f/F 


6.37 


Worst Sf/f 

Mean 

Worst 


31-h 

106.4 

1307 


Table 5 


0.C5 

1 

0.10 

7.21 

9*17 

79.9 

135.2 

83.6 

79.7 

13C^ 

1050 












